It is obvious that new physics is conducted now and willcontinue in rare processes in 
a foreseen future. Prior to their direct detection on colliders, supersymmetry effects may 
manifest themselves implicitly in rare processes. Supersymmetry gives elegant solutions 
to the problems of the standard model and predicts (minimal supersymmetry) the masses 
of supersymmetry particles in the range lOOGev — 0(l)TeV (with the exception of the 
lightest supersymmetry particle). In the present paper we consider the S°-meson's rare 
double-photon radiative decays B° — > 77 in the framework of the minimal supersymmetric 
standard model(MSSM). The process B s — > 77 is similar to the decay K° — > 77. However 
in the former process long-distance effects does not play an important role and short- 
distance contributions are the leading one. It should be mentioned that the two photons 
in the final state can be both in CP-odd and CP-even configuration. This possibility could 
become as a nontraditional source for the search for CP-violating effects in B-physics. 

The investigation of the i?-meson's rare decays is of great interest in order to test 
the standard model and beyond standard model physics. It seems that experimental 
investigation of beauty physics will be one of the major topics at the world facilities. The 
L3 collaboration(CERN) established that[l] 

Br{B s -> 77) < 14,8 • 1(T 5 

Br[B d -> 77) < 3, 95 • 1(T 5 (1) 

The experimental interest stimulates theoretical investigations [2-5]. All group of the 
authors have obtained (in the frame of the standard model) the same values for the 
branching ratios in the leading order 1/Myy : 

Br{B s -> 77) = (3.0 ± 1.0) • 10~ 7 

BrOB, ^77) = (1.2 ± 1.1) -10- 8 (2) 

In the framework of supersymmetric extension of the standard model one has following 
set of diagrams which contribute into decays B — > 77 (the diagrams are classified by 
particles which appear in the loops): a) charged gauge fermions(chargino) and up scalar 
quarks, b) charged Higgs particles and up quarks, c) neutral gauge fermions(neutralino) 
and down scalar quarks, d) gluino and down scalar quarks. The Lagranjian of interaction 
down quarks witth supersymmetric particles has the following form[6]: 
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L d - gd = -V~2g s ^ l -^dfg a d L3 + ^g^F^^dfgd^ (3) 

where (7 = esin^vK, V and [/ are charged gauge fermions' mixing matrices, C is charge 
conjugtion matrix, T L , T R are up skalar quarks mixing matrices, tg(3 = V2/V1, i>i,2 are are 
scalar Higgs fields' vacuum expectation [6], O is Kabibo-Kobayashi-Maskava matrix, N 
is neutral gauge fermins mixing matrix [6], Qd = —1/3, F L ,F R are down scalar quarks 
mixing matrices, g s is strong interaction coupling constant, A a are Gell-Mann matrices. 

One can write down the amplitude for the decays B — > 77 in the following form, which 
is correct after gauge fixing for final photons: 

T(B - 77) = eKhyMiAg^ + iBe^k^] (4) 

where e%(ki) and e^fo) are the polarization vectors of final photons with momenta k\ 
and k 2 respectively. 

Let us mention that our calculations are performed in the frame of Feynman-t'Hooft 
gauge and we use dimension regularization technique for divergent Feynmann diagrams. 
Only one particle reducible diagrams contain divergent parts. The divergent parts mutu- 
ally cancel in the sum of amplitudes and due to the GIM mechanism [7]. 

Ccontribution of the diagrams with charged gauge fermions into amplitudes A and B 
have the following forms: 

„ |2 m,m t \U j2 \ 2 Mj, 
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where 

m 2 (u L i) j, 5 — 12x + 7x 2 + (Ax — 6a; 2 )£na; 

x ' = M^*y /l(;r) = 2(1 - x y ' 

, /N 29 - 96a: + 111a; 2 - 44a; 3 + 6a;(4 - 9a: + 6a; 2 )£na; 
/2(X) = 6(1^ • (6) 

The higgs particles contribution into A and B have the following forms(contributions 
of other supersymmetric particles are significant small): 

A = i^^G^sieQnfM^OuOUfiixi) + cot 2 Pf 2 ( Xi ) - tan 2 P^j^ jhfa)}, 

B = ^^^G F f B (eQ D rY,O b ^O:Afl(x ^ ) + cot 2 (3f 2 (x ^ )+t^ (7) 

where fs is the 5-meson decay constant, Oj, being elements of Cabbibo-Kobayashi- 
Maskawa matrix, I = s, d, Qd = — 1/3 is charge of down quarks, Mb s ~ m& + m s , tan /3 = 
V2/V1, (vi t 2 being vacuum expectation of the Higgs bosons) [6] and we have introduced 
following notation: 

m?(ui) — 3x + 8x 2 — 5a; 3 + (Qx 2 — Ax)£nx 

Xl = W(H)' h[X}= 2(1 -a;) 3 ' 

„ , . 31a: - 84a; 2 + 69a; 3 - 16a; 4 + 6a;(x 2 - Qx + A)£nx 

M*) = ' (8) 

One can obtain contributions other supersymmetric particles from formulae (5) by fol- 
lowing way: 

Charged gauge fermions and right handed up scalar quarks: 



V n -> 0, 



U j2 | 2 m 2 (u R i) I V j2 



cos 2 (3 m b m s sin 2 (3 ' 
m(u Ll )^m(u m ), r£-T* (9) 
M(u Ll ) - M(H), M(x*) - m(«,), (10) 
Neutral gauge fermions and left handed down scalar quarks: 

V jX -> 2[sin w Nji ^— (i + Q d sin 2 w )^- 2 ] = ^, 

COS (7^ z 

Uj2 - -^3, Q« - 0, M(x*) - M(x°), 

m(^)^m(J Li ), r£-F*, (11) 
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One can obtain expressions for the diagrams with neutral gauge fermions and righ handed 
down scalar quarks from the contribution of neutral gauge fermions and left handed down 
scalar quarks by following way: 



m s m b 
Ll v/2.\/„ cos >' Lj V2M W cos f3 J3 ' 



'"' -N* -> -V2(sm9 w Q d N* - q/^^N*) ee t] 



y2M\Y cos p cosCV 

""' -Njs^tRj, m(d Ll )^m(d m ), r£-F,f, (12) 



v^M^ cos /3 

One can obtain expressions for the diagrams with gluino from the contributions of the 
neutral gauge fermions by following way: 



iV; 3 -0, <7-^,, ^-->/2, $-0. (13) 

Using equation (4) we find that the branching ratio of the decays B — > 77 can be 
represented as 

Br < B ^ = 32^[ 4 l^ 2+ 5 M| l B ' 2 l- (14 » 

With the aid of equation (4) we can find the relative fraction of CP-odd state of two final 
photons is given by 

4A 2 

5 = 4A 2 + (1/2)M|5 2 ^ 
Using formulae (4)-(15) one can estimate the branching ratio and parameter S of the 
decay B s — > 77, We have used the following set of parameters: Y tot {B s ) = 5 ■ 10~ 4 eV, 
m t = HAGeV, m b = 4, 8GeU, m s = 0, 5GeV, f Bs = 200MeB [8]. For small value of 
the Higgs particles masses (M(H) lOOGeU) and large tan/3 the super symmetric con- 
tribution into branching ratio of the deacay B s — > 77 are significant large than standard 
model estimate (Br susy (B s — > 77) ~ 10~ 6 > Br SM (B s — > 77) ~ 10~ 7 ). In the wide 
range of supersymmetric parameters (1 < tanf3 < 50, lOOGeU < M(H) < 350GeU) 
the supersymmetric contribution on the same level or more as standard model estimate 
(Br SM (B s -> 77) ~ 10- 7 ). 

The authors express their deep gratitude to N. Amaglobeli, T. Kopaleishvili, A. Tavkhe- 
lidze, Z. Garuchava, T. Sakhelashili and M. Xaindrava for support and interesting dis- 
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